Does the Two-Dimensional t-J Model have Hole Pockets? 
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We have calculated the high temperature series for the momentum distribution function rik of the 
2D t-J model to twelfth order in inverse temperature. By extrapolating the series to T — 0.2J we 
investigate the possibility of hole pockets in the t-J model. We find no indication of hole pockets at 
an electron density of n = 0.9 with J/t = 0.5 or J/t = 1.0. 
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The locus in momentum space of low energy single par- 
ticle excitations is of prime importance for what it can tell 
us about the nature of the low energy degrees of freedom 
of the two dimensional t-J model. Two alternatives have 
been widely discussed: a) a large Fermi surface consistent 
with Luttinger's theorem for a nearly half filled band and 
b) small hole pockets centered around (tt/2, 7r/2), similar 
to a lightly doped valence band in a semiconductor [Q. 
From the available experimental data which of these is 
correct for high temperature superconductors, if either, 
is still an open question. Angle resolved photoemission 
experiments and neutron scattering experiments gener- 
ally support a large Fermi surface. However, transport 
measurements are more consistent with a small number 
of positive charge carriers, in line with b) above. The 
motivation for hole pockets centered around (7r/2,7r/2) 
came from theoretical studies of a single hole in a Neel 
background More recently, hole pockets have been 
discussed theoretically |^|| and experimentally Q in re- 
gard to shadow bands. 

A detailed investigation of hole pockets requires the 
complete single particle spectral function. However, 
Duffy and Moreo j| have discussed the signature of hole 
pockets in nk for a spin density wave mean-field approx- 
imation to the Hubbard model Q . Within this approach 
a deep notch centered on (7r/2, 7t/2) develops in nk for 
low enough temperatures. Below we investigate the pos- 
sibility of this type of hole pocket in the 2D t-J model. 

To investigate the nature of the occupied states for the 
2D t-J model we calculated the high temperature series 
for the momentum distribution rik to twelfth order in 
inverse temperature. This is an extension of an earlier 
series calculation by Singh and Glenister H to eighth 
order in inverse temperature. Our series coefficients agree 
with theirs through eighth order. The Hamiltonian for 
the t-J model is given by 
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with the constraint of no double occupancy. The defini- 
tion of the single spin momentum distribution function 
is 



with n r — (cg CT c rcr ). The series coefficients are gener- 
ated in coordinate space with the Fourier transform to 
momentum space done exactly. The resulting series co- 
efficients in momentum space are exact functions of the 
wave vector. By extrapolating the series to low tempera- 
tures and using the exact momentum dependence of the 
series we can probe the momentum dependence of nk 
with high precision. 

To reach low temperatures we need to analytically con- 
tinue the series for n^. A standard way to do this is to 
use Pade approximants. For nk the straightforward ap- 
plication of Pades does not work very well. One way 
to improve the convergence of Pade approximants is to 
change the expansion variable to move the poles of the 
Pades farther from the origin || . Exactly what change of 
variable to choose is difficult to know for unknown func- 
tions. However, if the Pades converge well for a certain 
choice of parameters for rik, we can calculate the ratio of 
nk for this choice of parameters to nk for another choice 
of parameters, generally close by. This gives a much sim- 
pler function to extrapolate with the Pades and much 
better convergence. Once the extrapolation is done the 
desired result can be obtained by multiplying the ratio 
by the known function. For nk it's easiest to fix n, J and 
t and vary k. This is the method of extrapolation used 
to obtain the data below. 

Hole pockets are expected to be most prominent for the 
underdoped cuprates Q, S ^ 15%. To check this we fix 
the electron density as n = 0.9 (5 — 10%). The coupling 
constant ratio is less well known, but is typically chosen 
to be around J/t ~ 0.4. As a representative value we use 
J/t = 0.5. Since strong antiferromagnetic correlations 
are expected to promote hole pocket formation we have 
also investigated nk for J/t = 1.0. The results of our 
calculations are shown in Fig. 1. The momentum distri- 
butions for the two values of J/t are quite similar, with 
nk slightly larger near (tt, tt) for J/t — 1.0. Since for both 
values of J/t nk must satisfy the sum rule ^ k nk = n/2 
there must be a corresponding redistribution of weight 
for nk at J/t = 1.0 elsewhere in the Brillouin zone. To 
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check this would require a more complete set of data 
than presented here. The main feature of both curves 
presented in Fig. 1 is that rik for the parameters chosen 
here is a smooth function of k. Our results also satisfy 
the inequality S rik < (1 + <5)/2, with rik = 0.55 for 
k = and rik monotonically falling with increasing k . 

From the smooth, monotonic k dependence of rik there 
is clearly no indication of hole pockets near (ir/2,ir/2). 
There is also no obvious sign at T — 0.2 J of a quasi- 
particle discontinuity. Where to look for a quasiparticle 
discontinuity is a subtle problem when T ^ ||. From 
Fig. 1 we can see that rik = 1/2 is at a different wave 
vector than where rik has its maximum gradient, so we do 
not obtain a unique k F . This problem will be addressed 
in a future publication ||. 
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FIG. 1. The momentum distribution plotted versus 
wave vector k = (x7r,x7r) along the diagonal of the square 
Brillouin zone. The electron density is fixed at n = 0.9 and 
the extrapolated temperature is T/J = 0.2 (room tempera- 
ture), a) J/t = 0.5 b) J/t = 1.0. 

The momentum distribution for strongly correlated 2D 
electrons has been studied previously by a range of nu- 
merical methods. Stephan and Horsch || studied the 2D 
t-J model using the Lanczos algorithm at T = on 16-20 
site clusters. They had relatively few k-points available 
due to the small cluster sizes, but their results are clearly 
consistent with a large Fermi surface. The previous high 
temperature series work [pi presented results for a range 
of dopings at T — 1.0 J, all of which are consistent with 



a large Fermi surface. Quantum Monte Carlo has been 
used to study rik in the Hubbard model |l(]] . Results for 
n = 0.87, U/t = 4, and fit — 6 on a 16 x 16 lattice are 
also consistent with a large Fermi surface. 

Recently, the question of hole pockets in the Hubbard 
model has been reexamined |Q , with the suggestion that 
hole pockets would not show up until fairly low temper- 
atures of order fit ~ 10 for U/t ^ 8. While the t-J 
model and the Hubbard model are only directly compa- 
rable for J/t <C 1, our data at J/t = 0.5 can be approxi- 
mately compared to the Hubbard data using the relation 
J/t = 4t/U. For J/t = 0.5 we expect results similar to 
those for the Hubbard model at U/t w 8. The temper- 
ature of our results at J/t = 0.5 is T = Q.2J = Q.lt. 
This is low enough that we would expect to see an effect 
due to hole pockets of the type discussed for the Hubbard 
model (3) . Our current results do not show any indication 
of hole pockets of this type. 

In summary, we calculated a twelfth order high tem- 
perature series for the momentum distribution rik of the 
2D t-J model. By extrapolating the series to T = 0.2 J we 
examined rik near k = (7r/2,7r/2), looking for hole pock- 
ets centered at this wave vector. We find a smooth mono- 
tonic wave vector dependence for rik through (7r/2, 7t/2), 
with no indication of hole pockets. 
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